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Abstract 

In a previous paper we have proven that any multi-resolution analysis of £ 2 {R) 
produces, for even values of the inverse filling factor and for a square lattice, a 
single-electron wave function of the lowest Landau level (LLL) which, together with 
its (magnetic) translated, gives rise to an orthonormal set in the LLL. We have also 
discussed the inverse construction. 

In this paper we simplify the procedure, clarifying the role of the kq-representation. 
Moreover, we extend our previous results to the more physically relevant case of a 
triangular lattice and to odd values of the inverse filling factor. We also comment 
on other possible shapes of the lattice as well as on the extension to other Landau 
levels. 

Finally, just as a first application of our technique, we compute (an approxima- 
tion of) the Coulomb energy for the Haar wavefunction, for a filling v = ^. 
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I Introduction 



Since its discovery, the fractional quantum Hall effect (FQHE) has been considered as 
a very exciting problem to be understood. For this reason, it is not surprising that the 
presence of the plateaux in the Hall resistivity has been explained in many different ways 
by many different people. Most of these proposals are analyzed in monographies like 
[20, 6, 10]. The current belief is that the Laughlin wave function \p£ can explain the 
existence of these plateaux. As it is known, corresponds to an incompressible fluid 
and it describes electrons which are strongly correlated. It is also well know, however, 
that for certain values of the filling factor u, the electrons behave like a two-dimensional 
crystal, so that they must be described by a sustantially different wave function, 
which is usually taken as a normalized Slater determinant constructed starting with single 
electron normalized wave functions which should be mutually orthogonal in order vj/Cw) 
to be normalized in the thermodynamical limit, [3]. Many proposals in this direction has 
been made during these years, [22, 16, 18, 23, 14, 9] among the others, most of which 
are unable to reproduce the critical value of the filling factor for which a transition from 
liquid to crystal has been observed, [17, 11, 26]. 

In this paper we are mainly interested to those values of v for which the electrons 
behave as a 2-dimensional crystal. This aspect of the FQHE has already been considered 
in many papers, [22, 18, 14] etc.. What is new here is the use of a technique coming from 
wavelet theory, the multi- resolution analysis (MRA) of C 2 (R), which appears to be rather 
promising, as we will show in this work. 

In a previous paper, [2], we have proven that for a two-dimensional electron gas 
(2DEG) in a strong orthogonal magnetic field there exists a deep relation between any 
MRA of C 2 (R) and a normalized wave function for the N electrons system in the lowest 
Landau level (LLL). In particular we have seen that any MRA produces a single electron 
wave function which, together with its (magnetic) translated, produces an orthonormal 
set in the LLL. We have also proved that the procedure can be inverted. This surpris- 
ing result, however, was obtained assuming that the electrons live into a square lattice 
and, more important, was shown to hold only for a filling factor v = L G N. These 
mathematical constraints are physically rather unpleasant. In fact, classical and quan- 
tum arguments suggest that for low density the electrons should arrange themselves into a 
triangular lattice, [5, 3], rather than into a square one. This is due to the fact that the tri- 
angular lattice minimizes, for a fixed electron density, the Coulomb energy. Therefore, in 
order to have a deeper understanding about the concrete relevance of wavelets in FQHE, 
it is absolutely necessary to extend the results in [2] to the triangular lattice. Moreover, 
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the analysis of this new geometry, gives a new insight about the procedure and, in partic- 
ular, clarifies the role of the /cg-representation, which was an essential mathematical tool 
in our original approach. 

This is not the end of the story. It is well known that the main plateaux in the 
experimental measures of the resistivity tensor appear for filling factors like v — |, |, |, ... 
This fact seems to suggest that our technique, [2], is useless in most physically relevant 
situations (even if fillings like l/2n have a physical interest by their own right!). 

In this paper we overcome both these limitations: we show first how to obtain an 
orthonormality condition for the single electron wave functions for the triangular lattice. 
We also show that the same procedure can be adapted to lattices of arbitrary shapes, 
provided that the so-called rationality condition is satisfied. Incidentally, this extension 
reduces the relevance of the /cg-representation with respect to [2]. 

We also prove another result which looks quite fascinating to us: the odd values of 
z/ -1 are recovered by MRA of a different kind, that is MRA with dilation parameter d 
different from 2. In particular we will prove that a <i-MRA, d natural and bigger than 
2, produces an orthonormal set of single electron wave functions in the LLL related to a 
filling v = 2- The extension to other Landau levels is also discussed in some details, as 
well as the possibility of reversing the construction. 

The paper ends with an example of a 3-MRA which extends the usual Haar con- 
struction, [7], which is used to produce an iV-electrons Slater determinant made up with 
orthonormal single electron wave functions, for v — |. Finally, we compute the quantum 
Coulomb energy associated to this state under some simplifying approximations. 

II Mathematical tools 

In order to keep the paper self-contained we now quickly review, for reader's convenience, 
the main definitions of the mathematical tools we will use in the rest of the paper. 

II. 1 d-Multi-Resolution Analysis 

The main result in the theory of <i-MRA is the recipe which allows to construct an 
orthonormal basis in C 2 (R) starting from a set of functions i = 1,2, ...,d — 1, and 
acting on these with dilation and translation operators. 

The full story for d = 2 may be found, for instance, in [7], which also contains some 
remarks concerning d = 3. For a more complete reading on this subject we suggest 
references [21, 13], where the definition of MRA is extended to C 2 (R n ) and the dilations 
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and translations are not necessarily related to the integer numbers. Here, however, we 
give the definition of a <i-MRA for C 2 (R) and d integer (and d > 2), and then we focus 
only on those aspects which will be useful in the following. 

A d-multi-resolution analysis of L 2 (R) is an increasing sequence of closed subspaces 

... C V_ 2 C V_i C V C V 1 C V 2 C (2.1) 

with \Jj £Z Vj dense in C 2 (R) and HjezVj — {0}> an d sucn that 

(1) f{x) e Vj & f(dx) e V j+1 

(2) There exists a function <p G V , called a scaling function, such that {0(x — /c), A- G Z} 
is an o.n. basis of Vq. 

Combining (1) and (2), one gets an o.n. basis of V}, namely {(j)j :k (x) = d j / 2 (f)(d j x — fc), 
fc G Z}. 

For an integer d > 2 the existence of a set of functions ^W, % — 1,2, — 1, such 
that {V^fc^) = d? I 2 tl)( l \d? x — k),j,k £ Z i = 1, 2, d — 1} constitutes an o.n. basis of 
C 2 (R) can be proved. We are not interested in all the details of this construction here, 
see references [21, 13], but only to some easy consequences of this definition. 

Notice that the inclusion Vq C V\ yields the relation 

oo „ 

<f>(x) = \fd h n (f)(dx — n), h n = (0i, n |0) = Vd / <p(dx — n)(f)(x)dx. (2.2) 

n=-oo J R 

The orthonormality requirement of the functions 0, f R (f)(x — n)(f)(x)dx = 5 Ut Q, together 
with the equation (2.2), produces the key equation for the coefficients h n : 

^ hnhn+dl = $1,0, (2.3) 

which generalizes equation (2.12) of [2] to a dilation parameter d > 2. 

It is also easily checked that, if f R <p(x)dx ^ 0, then the coefficients must also satisfy 
the sum rule 

J2h n = Vd. (2.4) 

nez 

Equation (2.3) is all what we will need in the following. It is clear from the above 
construction that any <i-MRA produces a set of complex quantities h n which satisfies this 
condition. The converse procedure, which is not so interesting for our aims, is discussed 
in [21] and requires some care: in particular, it is easy to produce examples, for d > 3, 
of {h n } satisfying condition (2.3) and a scaling function <p satisfying the rf-scale relation 
(2.2), which do not produce a rf-MRA. 
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II. 2 kq-representation 

We give here only few definitions, namely those relevant for this paper, and refer to 
[24, 15, 25] for further reading and for applications. 

The genesis of the kq-representation lays in the well known possibility of a simultaneous 
diagonalization of any two commuting operators. In particular, the following distributions 

lO-jr 9-7T 

ipk q (x) = J— £ e ihna 6(x - q - na), k G [0, a[, ? e[0-[ (2.5) 
V a ngZ a 

are (generalized) eigenstates of both T(a) = e ipa and r(^) = e lx2lT l a . Here a is a given 
positive real number. 

These ^pkq(x) are Bloch-like functions corresponding to infinitely localized Wannier 
functions, [25]. They also satisfy orthogonality and closure properties. This implies that, 
roughly speaking, they can be used to define a new representation of the wave functions 
by means of the integral transform Z : C 2 (R) — > C 2 (V), V = [0, a[x [0, ?f[, defined as 
follows: 

h(k,q) := (ZH)(k,q) := / dulf^J) H (to) , (2.6) 

for any given function H(u) G C 2 (R). The result is a function h(k,q) G C 2 (V). 

To be more rigorous, Z should be defined first on the functions of C™(R) and then 
extended to C 2 (R) using the continuity of the map, [8]. 

Replacing ipkq(x) with its explicit expression, formula (2.6) produces 




h(k, q) = (ZH)(k, q) = \ — £ e~ lkna H (q + na), (2.7) 



a nez 



which can be inverted and gives the x-representation H(uj) G C 2 (R) of a function h(k, q) G 
C 2 (V) as follows: 

H(uj) = (Z' 1 h)(uj)= [ dkdqip kq (u)h(k,q). (2.8) 

Jv 

Due to (2.5), this equation can be written as 

h-ir r b 

H{u + na) = J— / dke ikna h(k,u), Vo; G [0, a[, Vn G Z. (2.9) 
| a Jo 

III Stating the problem 

The physical system we are considering is a 2DEG, that is a gas of electrons constrained in 
a two-dimensional layer belonging to the (O; x, y) plane, in a neutralizing positive uniform 
background and subjected to an uniform magnetic field along z. 
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The hamiltonian of the iV-electrons system can be written as 

HW = H™ + \(HW + H™) (3.1) 
where is the sum of iV contributions: 

^r = E^o(0- (3-2) 

1=1 

Here H (i) describes the minimal coupling of the electrons with the magnetic field: 

//o = i( £ + 4(r)) 2 = i( P ,-|) 2 + i( ft + |) 2 . (3.3) 

Notice that we are adopting here the symmetric gauge A = |(— y, x, 0) and the same units 
as in [3]. is the canonical Coulomb interaction between charged particles: 

and is the interaction of the charges with the background, whose explicit form is 

irrelevant here and can be found in [3]. 

In the following we will consider, as it is usually done in the literature, \(H^ + 
H^) as a perturbation of the free hamiltonian Hq N \ and we will look for eigenstates 
of Hq N ^ in the form of Slater determinants built up with single electron wave functions. 
This approach is known to give good results for low electron (or hole) densities, [18, 3]. 
However, comparison with the experiments shows also that the wave-functions proposed 
in [3, 18, 22, 9] does not reproduce completely the experimental data, as far as the 
transition from the Wigner crystal to the Laughlin liquid is concerned. Other wave- 
functions proposed in the literature share with this the same problem, so that the problem 
of finding the correct single-electron wave function describing the system for low electron 
density is still open. In this paper we will show how an o.n. set of wave functions in 
the LLL (as well as in any other Landau level) can be constructed easily starting from a 
d-MRA. 

The easiest way to attach this problem consists in introducing the new variables 

P' = Px -y/2, Q'= Py + x/2. (3.5) 
In terms of P' and Q' the single electron hamiltonian, H , can be written as 

H,= l -(Q' 2 + P' 2 ). (3.6) 
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The transformation (3.5) can be seen as a part of a canonical map U from (x,y,p x ,p y ) 
into (Q, P, Q', P') where 

P = p y - x/2 + ^={p x + y/2), Q = Px + y/2. (3.7) 
These operators satisfy the following commutation relations: 

[Q, P] = [Q', P'] = i, [Q, P'} = [Q', P] = [Q, Q'] = [P, P'] = 0. (3.8) 

Using [19], we deduced in [4] that, as a consequence of the above canonical transformation, 
a wave function in the (x, ?/)-space is related to its expression in terms of the new variables 
(P, P') by the formula 

iy / V \ 

iKs, y) = r r e^-^ P + PP '-^m P') dP dP>. (3.9) 

We want to stress that this formula differs from the analogous one in [2] since the map U : 
(x, y, p X) p y ) — > (Q,P,Q* ', P') differs from the one used in [2] . The reason for this difference 
is in the shape of the lattice, which produces a different canonical transformation. This will 
appear clearly in the following, when we will discuss how the lattice should be constructed 
explicitly and the role of the magnetic translations in this construction. 

The usefulness of the new variables stems from the expression (3.6) of H . Indeed, in 
this representation, the single electron Schrodinger equation admits eigenvectors ip(P, P') 
of H of the form i/>(P,P') = f(P')h(P). Thus the ground state of (3.6) must have the 
form f (P')h(P), where 

/o(P') = ir-V'e- 1 "' 2 , (3.10) 

and the function h(P) is arbitrary in C 2 (R), which manifests the degeneracy of the LLL. 
With /o as above formula (3.9) becomes 

il>(x,y)= r K®(r,s)h(s)ds, (3.11) 

J — oo 

where is what we call the kernel of the transformation, 

K^ir, s) = e \ * jvs-Pix-k+s? (312) 

where (3 = |(1 — ^). Here we use the suffix t to emphasize the shape of the lattice, 
triangular in this paper, while the index means that we are working in the LLL. Just 
to compare this result with the one obtained for the square lattice, [2], we recall that 
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Kq\l,s) = ^^3/4 e iys e ( x+s ) 2 / 2 . The extensions to higher Landau levels (see again [2] 
for K[ s) ), can be found simply replacing f in (3.10) with the excited harmonic oscillator 
eigenstates: using f\ will produce, clearly, a wave function in the 1-th Landau level. 
However, for large magnetic fields, it is well known that only the very first levels are 
relevant in the analysis of the FQHE, and this is the main reason why, quite often, only 
the LLL is considered. This is also our choice here, and for this reason, along this paper, 
we will consider only Kq \ but for few comments. 

Once the generic form of any function in the LLL has been given, formula (3.11), it 
remains to choose properly one of them, that is, to fix h(s) by requiring that it minimizes 
the Coulomb interaction. To achieve this aim we will first construct a particular ground 
state of the (infinitely degenerate) free iV-electrons hamiltonian Hq N \ We use a suggestion 
coming from the classical counterpart of this quantum problem. It is very well known 
that the ground state for a classical 2DEG is a (triangular) Wigner crystal: the classical 
electrons are sharply localized on the sites of a lattice whose lattice spacing is fixed by 
the electron density, [5]. What we expect, and what was proven in [3], is that, at least for 
certain values of the filling factor, the quantum ground state should not be very different 
from this classical picture. 

Let us introduce the so-called magnetic translation operators T(ai) defined by 

T(al) = exp (iU c ■ a*), i = 1, 2, (3.13) 

where H C = (Q, P) and are the lattice basis vectors: 

ai = a(l,0), a 2 = ^(l,V3). (3.14) 

We assume the following (minimal) rationality condition on the area of the fundamental 
cell of the lattice: 

a lx a 2y - a ly a 2x = 2ir, (3.15) 
which for the triangular lattice reads 

a * = ^ (3-16) 

Notice that this is different from the analogous condition for a square lattice, a 2 = 2ir, 
[2]. Because of this condition we have 

[T(a-[),T(a 2 )} = 0, (3.17) 

which, together with 

[T(oi),# ] = [T(a 2 ),H ] = 0, (3.18) 



8 



which easily follow from (3.8), shows that T(dl),T(a 2 ) and H Q are a set of commuting 
operators. 

This implies that if tp is an eigenstate of H with eigenvalue e, then T(a 1 ) n T(a 2 ) m ^> is 
still an eigenstate of H corresponding to the same eigenvalue, for all n, /// G Z. This can 
be seen as another evidence of the infinite degeneracy of the different Landau levels. 

Because of (3.14), the magnetic translations take the form 

7\ := T(ai) = e <0 (P"+*) = e iaQ ', T 2 := T(a 2 ) = e *S((p-+«/2)+^-*/2)) = e <^P = 

(3.19) 

which show the relevance of the definitions in (3.7). Notice that the last equality for T 2 
above is a simple consequence of the rationality condition (3.16). 

Using this equality we can check that the action of 7\ and T 2 on a generic function 
f(x,y) G C 2 (R 2 ) is given by 

f m , n (x,y) := T™T?f(x,y) = (_i)™ e i(W-x„ mW)/(a . _ ^ y _ ^ (3>20) 

We see from this formula that if, for instance, f(x,y) is localized around the origin, 
then fm,n{x,y) is localized around the lattice site (X nm ,Y nm ). Here X nm = —a(n + y) 
and Y nm = —am^ = —m^-. This result extends, as expected, the one for the square 
lattice, [2], and explains why we call Tj translation operators. 

Moreover, equation (3.20) is the key relation by means the quantum triangular lattice 
can be constructed: we simply start from the single electron ground state of H given 
in (3.11), ip(x,y). Then we construct a set of copies ipm t n(x,y) of ip as in (3.20), with 
///. // G Z. All these functions still belong to the lowest Landau level for any choice of the 
function h(P), due to (3.18), and are localized around the sites of our triangular lattice. 
iV of these wave functions ip m ^ n {.x^y) are finally used to construct a Slater determinant 
for the finite system: 



^ (JV) (Ei,E2,-,EAr) 



1 

W\ 



■0mi,ni(ZLl) 
^m2,«2 (— l) 



1>- 



m N ,n N 



(hi) i> 



m N ,n N 



(L 2 ) 



'4 , m2,n2 {T-n) 



1>: 



m N ,n N 



[L N ) 



(3.21) 



We have already discussed in [2] the relevance of the ortonormality of differently localized 
single electron wave functions when the thermodynamical limit of the model is considered. 
For this reason we look for o.n. single electron wave-functions. Therefore, let vp(x,y) be 
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as in (3.11) and ip n>m (x,y) = T^T^^>{x,y). We are interested in finding conditions on 
h(P) such that the orthonormality condition (ONC) 

S n ,m-= ip ,o{x, y)ip n ,m{x,y)dxdy = 5 mfi 5 nfi , (3.22) 

J R 2 

is satisfied. Any function ip(x,y) which satisfies this condition and which belongs to the 
LLL can be used to generate the function ^ < - Ar - ) (r 1 , r 2 , ...,r N ) as in (3.21). Moreover, it is 
not hard to prove that the ONC can be rewritten as 

S n ,m ■= I dse ism ^h(s)h(s - na) = 5 mfl 5 nfl . (3.23) 

This is a consequence of eqs. (3.11), (3.12), (3.16) and (3.20). 

It is interesting to remark that this equation implies that ?/>o,o is normalized in £ 2 (R 2 ) 
if and only if h(s) is normalized in £ 2 (R). This is clearly a consequence of the canonicity 
of the map U. 

Moreover, equation (3.23) coincides exactly with the analogous one we have obtained 
for the square lattice, under the condition a 2 = 2tc, [2]. Incidentally, this fact strongly 
suggests that the same condition will appear independently of the shape of the lattice. 
Finally, we observe that the ONC (3.22) produces the same equation (3.23) independently 
of the particular Landau level we are considering. In other words, it is not hard to check 
that we obtain equation (3.23) even if the kernel K<p is replaced by Kf\ for any given 
I E N. This proves, therefore, that the problem of the orthonormality of the single 
electron wave functions is Landau-level independent. Of course, this does not means that 
the o.n. wavefunctions look the same in any Landau level: this is true only at the level 
of the function h(s). On the contrary, the wavefunctions in the coordinate representation 
will be different because the kernels of the transformations, K^\r,s), are different for 
different values of l\ 

If we require equation (3.23) to hold for all integers n,m we are implicitly assuming 
that all the lattice sites are occupied by an electron, so that the filling is v — 1. 

In [2] we have discussed how the ONC should be rewritten in order to cover a fractional 
value of the filling factor. A possible way for requiring orthonormality between single 
electron wave functions for v = |, d G N, consists in replacing (3.23) with the following 
equation, 

/ dse isdm ^h(sjh(s - na) = 5 mfi 5 n , , (3.24) 

for all ///. // G Z, which is obtained under the assumption that only one lattice site every d 
is occupied. We see, therefore, that the ONC becomes very much filling- dependent, and, 
in a certain sense, this fact is not very surprising to us: a filling v < 1 can be interpreted as 
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if in the physical system we have less electrons than for v — 1, so that these electrons are 
more distant from one another in order to minimize the Coulomb energy. For this reason 
we expect different analytic expressions for the single electron wave functions depending 
on the value of the filling. 

We are left now with the main problem: is it possible to produce solutions of equation 
(3.24)? This problem has been partially solved in [2] where solutions where explicitly 
built up for even values of the inverse filling factor starting from a 2-MRA and making 
use of the /cg-representation. In the next section we will show how to solve (3.24) both 
for odd and even values of u^ 1 , and this will be done without even mentioning the kq- 
representation. However, in Section V, in order to show the relation with our previous 
approach, we will rewrite everything using the Zak transform. 



IV d-MRA and odd inverse filling 

In this section we propose a method for constructing solutions of equation (3.24) which 
works for all integer values of the inverse filling. This method produces exactly the same 
solutions as in [2] only for v — |, while produces different wave functions for other even 
values of v~ l . Moreover, it will produce solutions also for those values of v~ x which where 
not covered by the approach in [2], that is for odd values of v~ x . 

Let us consider a given ci-MRA of C 2 (R) and its related square-summable set of com- 
plex numbers {/i n }nez satisfying condition (2.3). Following [2], we use these coefficients 
to define a function T d {u) as follows 

f ^tY.igz hie ils ^, s G [0, a[ , x 

T d (s) = { ^^ lez ' L ' L (4.1) 

j 0, otherwise. 

It is clear that Td(s) is square integrable and not periodic. In particular, due to condition 
(2.3), we have ||T d ||l = J R \T d (s)\ 2 ds = 1. 

With this definition it is straightforward to check that 

/ dse lsdm ^Us)T d (s - na) = 5 n0 5 m0 . (4.2) 

JR 

Infact: 

a) the supports of T d (s) and T d (s — na) are disjoint but if n = 0; 
b) 

/ dse^ dm ^\T d (s)\ 2 = - hh l a dse^ dm - l+ ^ = £ h p h p+dm = 6 m0 , 
jR a i, P ez Jo P ez 
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due to condition (2.3). 

The physical consequences of this fact are evident: calling 



i>d(x,y)= Kq\l, s)T d (s)ds 

JR 



(4.3) 



then the following ONC 



<^,T"Tf> d >=5 n0 S m0 , 



Vn, in 6 Z, 



(4.4) 



holds. We must stress that our construction does not imply that ipd is orthogonal to 
T^T^ipd for all integer values of i and j, but only for those indices which belong to 
the set { (Z, dZ) } . 

The set {T^T^tpd^jV) '■ n,m G Z} generates a triangular lattice where not all the 
sites are occupied: with our choice, only a site every d along the a 2 direction is non 
empty while all the sites along d[ are full. It is evident that this is only one among 
many choices: another possibility is exactly the opposite one: all the sites are occupied 
along 02 while only one every d is non empty along d[. This choice corresponds to the 
set {T^T^t/jdixiy) : n,m e Z\. More symmetrical choices are not difficult to imagine. 
Clearly, from a mathematical point of view, all these choices are equivalent. However, 
physical reasons suggest that numerical differences could arise in the computation of the 
Coulomb energy, the reason being the existence of possible different asymptotic behaviors 
of the single electron wave functions for different choices of the above sub-lattices. An 'a 
priori' decision seems very hard: the only reasonable criteria we could use are suggested 
by symmetry considerations. However we will see in Section VI that already the choice 
discussed above, even if it breaks down the symmetry for rotations of n/3 which one would 
expect for the triangular Wigner crystal, produces a fairly good localized eigenfunction 
of Hq and, therefore, is surely worth of a deeper analysis also form a numerical point of 
view. 

REMARKS.- (1) It can be useful to remark that this procedure is not unique, that 
is it produces as many o.n. sets in the 1-th Landau level (for any given shape of the lattice) 
as many d-MRA we are able to produce. This is not surprising because it reflects once 
more the infinite degeneracy of the Landau levels. Of course, among all the possibilities, 
we are interested in finding the one which minimizes the Coulomb energy. 

(2) We observe that our strategy suggests some kind of transition related to the value 
of v. The point is the following: in [3] our lattice was essentially constructed with a 
superposition of gaussian functions, localized around different lattice sites. The same 
single-electron wave function was to be used for any value of the filling factor, the only 
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difference being that the different wave functions were more or less distant depending on 
v, which, as in the classical case, appears essentially as a parameter fixing the distance 
between two neighbouring electrons. Here the situation is more delicate. Infact, even if 
the expression for T d is formally always the same independently of the filling, see (4.1), 
the explicit expression changes because the coefficients h n for a, say, 2-MRA are in general 
different from those of a 3-MRA or, in general, of any other <i-MRA. We will show this 
in Section VI, considering as an example the Haar 3-MRA. The consequence of this fact 
is that a <i-MRA produces an A^-electron wave function ip d N \r_\-,L2i ■■■■>Ln) an d two such 
functions, and for instance, need not to have the same analytical dependence on 

3 t 5 

their variables. A natural question to be considered is the following: is there any relation 
between a given il) d N \r_iiL2i ■■■■>Ln) an d the plateaux of the Hall resistivity corresponding 
to v — ^? We hope to be able to consider this question in a close future. 

As in [2], we show now just for completion how the procedure can be inverted, that is 
how a set of coefficients satisfying equation (2.3) can be recovered by a function satisfying 
the ONC (4.2). The idea is quite simple: given such a function K d (s) we simply define 




and it is now an easy computation to check that 

£ H n H n+dl = 5 l0 , VZ e Z. (4.6) 

neZ 

Indeed, using formula J2n^z emx ~ — a Hn&z&{ x ~ na ) we have: 




Our claim now follows from the assumptions on K d . 

Finally, it is not a big surprise that whenever we take K d (s) coincident with the 
function T d in (4.1), then the coefficients H n in (4.5) coincide with the h n , h n = H n for 
any n. 

Before ending this section, it may be worth stressing that the definitions (4.1) and 
(4.5), which are given here without any reasonable justification, are really a consequence 
of the analysis produced in [2], and, in particular, of the use of the /cg-representation, 
which was crucial to make evident the relation between MRA and ONC. 
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V The role of the /cg-representation 



In the previous section we have discussed the relation between a <i-MRA and an o.n. 
set of single electron eigestates of the free Hamiltonian of a 2DEG in a given Landau 
level. Our main requirements are that we want to generate a triangular lattice and that 
we want our approach to be useful for describing those values of the electron densities 
for which the plateaux are observed in experiments. These results, which extend in a 
relevant way the ones discussed in [2], have been obtained here without even mentioning 
the /cg-representation. Nevertheless, as we have just remarked, this has proved to be a 
crucial guideline in [2] . For these reasons, and because the ONC can again be written in a 
much simpler form in the variables (k, q), we devote this section to restate our results by 
making use of the Zak transform. Moreover, this approach will also explain why the role 
of the geometry of the lattice and the particular Landau level are completely irrelevant 
as far as we are only interested to the orthonormality between the vectors ip and T™T™tp. 

We begin with adapting the definitions sketched in Section II to our situation, that is 
considering the unitary operators Ti and T 2 defined in (3.19), Ti = e lQa and T 2 = e lP ^~ . 
A simple extension of the original Zak's proof produces the following result: 

Lemma 5.1 

The following set 

^^l/fE^Mf-Hna), ke[0,a[, q e [Q,—[ (5.1) 

satisfies 

T^ kq (p) = e^ kq (p), T 2 ^ kq {p) = e ik ^^ kq (p), (5.2) 

and 

(■a r2ir/a 

i) ^ JO ^^qiP^kqi^) = 5(p -p'). (5.3) 



This means that any ip kq (p) is a common eigenstate of T\ and T 2 for all k and q and, 
moreover, that they form a complete set in C 2 (V), V = [0, a[x[0, ^[, as in Section II. 

The completeness of this set can be used to write S njm in a different way. If we call 
h(k, q) the Zaq transform of h(s), 

h(k,q) = (Zh)(k,q) := ^^e^ + na), (5.4) 

V a nez 



we can find that 



S n , m = I dkdqe ma ^ km ^\h(k,q)\ 2 , (5.5) 
Jv 
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which coincides essentially with the one in [2]. Incidentally, once more we have an evi- 
dence of the canonicity of maps in the following sense: \\ip(x, y)\\c 2 (R 2 ) = HM s )ll£ 2 (R) = 
\\h(k, q)\\c 2 (v)- Requiring orthogonality only among ^oo and VVi,dm, n,m G Z, we can 
repeat essentially the same steps as in [2] and we conclude that S ntdm = ^no^mo if and 
only if 

Uk,q):=p\h(^,q)\* = ± (5.6) 

which looks more friendly than (3.24) because no integration appear! 

The solution of this equation can be obtained starting from a set of coefficients {h n } 
arising from a <i-MRA, and therefore satisfying (2.3). It is enough to define the function 

t d (k, g) = -L £ h n e ikn ^, (k, q) E V. (5.7) 
V 2tt neZ 

It is not hard to check now that (1) t d (k,q) satisfies condition (5.6), and (2) the inverse 
Zak transform of t d (k,q) returns T d as in (4.1). Indeed we have: 



J d (k,q) :=- V h.W l e«- t > 3 3}(l + j<- t > 3 f + + + e^^) = 

27r Mez 

s, P ez j=o 

Since Y^Zli^^V is equal to d whenever p — 0, ±d, ±2d, ±3d, ... and to for all the other 
values of p, we conclude that, with a change of variable m = dp, 

J d (k,q) := ^- h s h s+md e- im2 ^ , 



A. 

As for our second claim, we have 



which is finally equal to ^ as a simple consequence of equation (2.3). 



Z H d )(s)= / ip kq (p)t d (k 1 q)dkdq 



v 



fn __ ra /-27r/a 9 

— Y, h n dk dqe iqma 5(p — k + ma)e ink - 

-^J2 h n [ a dk5(p-k)e mk ^, 
\/a frt, Jo 



neZ 



which is equal to zero whenever p [0, a[ while gives a non trivial contribution otherwise. 
Computing the integral for p e [0, a[, our assertion easily follows. 
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In analogy with our results in the previous section and with those in [2] , it is not hard 
to imagine how to use a function ta(k, q) satisfying the ONC in its form (5.5) to generate 
a set of coefficients H n satisfying equation (2.3), reversing in this way the procedure and 
obtaining a <i-MRA from an o.n. set in the LLL. 

We want to conclude this section with a remark concerning the fact that the ONC 
looks identical for the square and for the triangular lattice and also for all the Landau 
levels different from the lowest one. These facts can be nicely understood in terms of the 
/cg-representation. The idea is the following: 

let ip be an eigenstate of the self-adjoint operator H belonging to a certain eigenvalue 
e, and let us suppose that H commutes with the two unitary operators T 1? T 2 defined as 
in (3.19). It is clear that, putting ^ nm = T™T™tp, then H fy nm = e\l/ nm , for all integers n 
and m. We want to answer now to the following question: 

is it possible to fix tp in such a way that 

S m ,n ■ = < 1p, ^nm >= $nO$mO (5.8) 

holds? 

The answer follows quite easily now from Lemma 5.1. In fact, using the completeness 
of the set ipkq associated to the magnetic translations, we can write 

S m ,n = I dkdq<i),i) kq ><i) kq) T™T 2 "V >, 
Jv 

and we go back to (5.5), putting h(k,q) =< ipkq,ip >> simply because each ip kq is an 
eigenstate of both Ti and T 2 . As we see, the role of the shape of the lattice, as well as the 
particular Landau level we want to consider, are completely unrelevant for this result as 
we have explicitly seen in Section III. 

VI An Example: Haar for v = \ 

In this section we will discuss the easiest non trivial example of a 'trial' ground state for 
the 2DEG arising from a 3-MRA with an Haar-like behavior. Let us first introduce the 
3-MRA. This can be done as follows: the characteristic function of the unit interval [0, 1[, 
H(x), is known to be a scaling function for a 2-MRA with coefficients ho — h± — 
hi = for all the other values of %. The same function can be easily shown to be a scaling 
function also of a 3-MRA with coefficients h = hi = h 2 = hi = for all the other 
values of i. It is also obvious that equation (2.3) is verified. Therefore, we use this set of 
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coefficients to construct our single electron wave function as shown before: the function 
T 3 (s) must be defined as 

[0, otherwise. 

Using T 3 into (3.11), together with the usual definition of the error function, [12], we 
obtain 

= ir o„ \ \ o mS g{ ?> y > y) + G( - x > y>y+—) + °( x ' y>y + — > (°) 

J6a(3(l + ^)2tt 1 /4 a a 

where, as before, we have defined /3 = |(1 — ^) and 

G(x,y,a) = e -fg- <a (*-»/^ [$( ^ ~ ~ m + 2a ^ ) _ ~ ^M) ~ 

(6.2) 

■03(e) is now used to define (V^nmCn) as in (3.20), and one can explicitly check that 



f{2 V'3(£)(V'3)n3m(£) d L = 5 n0 5 m0 . 

An interesting feature of these functions is their asymptotic behavior which can be found 
considering the asymptotic behavior of the error function, [12]. After some algebra we 
obtain: 

^ W v^vr 3 / 4 X 2P(x - y/y/3) -iy 

_| \ | \ e -l3a 2 -a(2f3(x~y/V3)-iy) 



2(3{x - y/y/3) - i(y + 2ix/a) 2(3{x - y/y/3) - i(y + 4ir/a) 

1 1 

+ —, — ^ ; r- + 



2(3(x - y/y/3,) -iy + 2a[3 2(3{x - y/y/3) - i(y + 2tx/o) + 2a(3 

+ t= II- (6.3) 

2(3{x - y/y/3) - i{y + An /a) + 2a(3 s 

What is interesting about this formula is that it shows that whenever x — y/y/3 ^ 
then 03 (r) decreases very fast while, if x — y/y/3 = 0, then the decay is rather slow (like 
1/| y | and independent of x). This is rather appealing because suggests that, even if along 
a direction the decay is the one expected by the Balian Low theorem, [1], outside this set 
of zero measure the wave function decreases very fast in both variables. 

Remark.- In [2] we have discussed an example for v — \, also related to the Haar 
basis. We give here the result of that analysis in order to compare old and new results 
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and, at the same time, to correct two typos. The wave function generating the o.n. set is 

whose asymptotic behavior can be found with the help of [12]: 

T 2 (x,y ~ V , 5/4 - + 

47r t '/ 4 x — xy 

_| 1 c -7r-o(x-iy) / 1 | ]_ \\ 

x — i(y — a) x + a — iy) x + a — i(y — a) 

In particular we see that the decay of this function is rather slow, but for a set of zero 
measure, when compared with that of ^3. 

Let us now use this wavefunction in the computation of the Coulomb energy. We 



<flW> 



repeat the same steps as in [3]. We want to compute the limit liniAr^oo — , where 

ip( N *) is constructed as in (3.21). We write this mean value as a sum of two terms, a kinetic 
contribution coming from H , which is very well known since all the single-electron wave 
functions are eigenstates of H , and a Coulomb correction E c which is the one we really 
need to compute. Therefore 

< HVO>,„ = + ^ Er = Ew + SE {6A) 

N~ >oo _/V 2 

where E w is the classical energy per electron of the Wigner crystal, E w = — 0.7821y / z/, 
[3], while SE is the quantum correction, 

6E=hml {[E d ((n,m),0)-E ex ((n,m),0)\-— —}. (6.5) 

N ^°° 2 (n,m)eC \^n,m\ 

Here C is the subset of Z x Z corresponding to the sublattice we are considering. In this 
case, for instance, C = (Z, 3Z)\(0, 0) to avoid self-energy divergences. Moreover, as in 
Section III, we have R n m = (X n>m , Y n>m ), while the direct and the exchange contributions 
are: 

|VVn(ri)lWr 2 )| 2 

and 



E d ((n,m),0) := d r 1 d r 2 . . , (6.6) 



E ex ((n,m),0) := / d r 1 d r 2 ; ; . (6.7) 

J |£i — £2! 
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In this paper we are not interested in giving a detailed result for SE, since the CPU time 
required to compute each integral is very high (more than 1 day for each integration with 
a reasonable approximation). For this reason, we only compute here the integrals using 
a Montecarlo method with 250000 points and restricting the integration to a large but 
compact region (taking advantage from the decay behavior of V\i,m an d checking first the 
normalization of each wavefunction). With this very rough approximation the CPU time 
for each integration is about four hours. Moreover, we disregard the contributions coming 
from the exchange integrals, which are usually some orders of magnitude smaller than the 
direct one, [3]. Finally, due to the fact that the direct energies Ed((n,m),0) go to zero 
when (n, m) increase, we also work with a finite (but big) lattice. 

For the above reasons we do not expect that the present result could have any real 
meaning, but nevertheless it is useful because gives a first idea of the numbers appearing 
in the game. The result we have obtained is SE = 0.3184 which must be compared with 
the much lower result in [3], 0.0657. It is useful to stress that this does not imply at 
all that the Haar wave-function is physically useless, but only that a different numerical 
integration technique has to be used. A more detailed numerical analysis of the procedure 
is now being performed. This will include the computation of the Coulomb energy for 
different values of the filling factor as well as the use of different d-MRA. We hope to be 
able to produce physically relevant results in a close future. 

VII Outcome 

In this paper we have carried on the analysis of the connection between a MRA of C 2 (R) 
and the FQHE, first proposed in [2]. In particular we have shown how a single electron 
wave function which, together with its magnetic translates, produces an o.n. set in the 
LLL can be constructed starting from a rf-MRA. This procedure works for v = |, d E N 
and d > 2. We have also shown that this procedure can be essentially inverted since to 
any o.n. basis of translated functions of the LLL (corresponding to v — |) corresponds a 
set of coefficients satisfying the main condition of a d-MRA of C 2 (R). 

All of these results have been obtained working with a triangular lattice, which is the 
one suggested by classical and quantum energetic considerations. The extension to higher 
Landau levels has also been briefly outlined. 

We have also given an example of our construction for v — |, using the Haar 3- 
MRA, and for the related o.n. basis in the configuration space we have also obtained an 
approximated value of the Coulomb energy. 

It is clear that our method produces a wide set of possible o.n. functions in the LLL, 
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and our future interest consists mainly in going beyond the simple Haar choice trying to 
minimize the Coulomb energy so to recover completely the experimental data. 

As already mentioned, it will be very interesting also to analyze the (possible) connec- 
tion between the different wave functions produced by our method and the plateux 
corresponding to v — |. Of course, this mechanism is completely non-standard, so that 
an analysis of the relations between our approach and the one originated by Laughlin 
should also be worked out. 

Acknowledgements 

I would like to thank Prof. M. Zarcone, Prof. F. Oliveri, Prof. G. Russo and Dr. F. 
Agostaro for their suggestions concerning the numerical computation in Section VII, and 
Prof. M. Sammartino for giving me most of his CPU time (for free!!!). 

I also would like to acknowledge financial support by the Murst, within the project 
Problemi Matematici Non Lineari di Propagazione e Stabilitd nei Modelli del Continuo, 
coordinated by Prof. T. Ruggeri. 



References 

[1] Antoine J. P. and Bagarello F., Localization properties and wavelet-like orthonormal 
bases for the lowest Landau level, in Advances in Gabor Analysis, H.G. Feichtinger, 
T. Strohmer Eds., Birkhauser, Boston, to appear 

[2] Bagarello F., Multi-Resolution Analysis and Fractional Quantum Hall Effect: an 
Equivalence Result, J. Math. Phys., 42, 5116-5129, (2001) 

[3] Bagarello F., Morchio G. and Strocchi F., Quantum Corrections to the Wigner Crys- 
tal. An Hartree-Fock Expansion, Phys. Rev. B 48 (1993) 5306 

[4] Bagarello F., Electrons in a magnetic field, Wigner Crystal and the Fractional Quan- 
tum Hall Effect, Ph.D. thesis SISSA/ISAS, year 1990-91 

[5] Bonsall L. and Maradudin A. A., Some static and dynamic properties of a two- 
dimensional Wigner crystal, Phys. Rev. B 15 (1977) 1959 



20 



[6] S. Das Sarma and A. Pinczuk, Perspective in Quantum Hall Effects, Wiley, New 
York, (1997). 



[7] Daubechies I., Ten Lectures on Wavelets, Society for Industrial and Applied Mathe- 
matics, Philadelphia, 1992 

[8] Daubechies I., Grossmann A., Meyer Y., Painless nonorthogonal expansions, J. Math. 
Phys., 27 (1986) 1271 

[9] Ferrari R., Fractional Quantum Hall Effect: construction of the Hartree-Fock state 
by using translational covariance, Int. Jour. Mod. Phys. B, 8, (1993), 529, Ferrari 
R., Two-dimensional electrons in a strong magnetic field: a basis for single-particle 
states, Phys. Rev. B, 42, 4598, (1990) 

[10] Girvin S.M., The Quantum Hall Effect: Novel Excitations and Broken Symmetries, 
Springer Verlag, Berlin and Les Editions de Physique, Les Ulis (2000). 

[11] Goldman V.J., Santos M., Shayegan M., Cunningham J.E., Evidence for two- 
dimensional quantum Wigner crystal, Phys. Rev. Lett., 65, (1990), 2189 

[12] Gradshteyn I.S. and Ryzhik I.M., Table of Integrals, Series and Products, Academic 
Press, New York and London 1980 

[13] Grochenig K., Madych W.R., Multiresolution Analysis, Haar Bases, and Self-Similar 
Tilings of R n , IEEE Trans. Inf. Th., 38, No.2, (1992) 556 

[14] Yang K., Haldane F.D.M., Rezayi E.H., Wigner crystal in the lowest Landau level at 
low filling factors, Phys. Rev. B, Rapid Commun, 64, 081301 (2001) 

[15] Janssen A.J.E.M., Bargmann transform, Zak transform, and Coherent states, J. 
Math. Phys., 23 (1982) 720 

[16] Lam P.K., Girvin S.M. Liquid-solid transition and the fractional quantum Hall effect, 
Phys. Rev. B 30 (1984), 473 

[17] Levesque D., Weis J.J., MacDonald A.H. Crystallization of the incompressible 
quantum-fluid state of a two-dimensional electron gas in a strong magnetic field, 
Phys. Rev. B 30 (1984), 1056 

[18] Maki K., Zotos X., Static and dynamic properties of a two-dimensional Wigner crystal 
in a strong magnetic field, Phys. Rev. B 28 (1983), 4349 



21 



[19] Moshinsky M. and Quesne C, Linear Canonical Transformations and Their Unitary 
Representations, J. Math. Phys. 12 (1971) 1772 

[20] R.E. Prange and S.M. Girvin eds., The Quantum Hall Effect, Springer Verlag, (1990). 

[21] Strichartz R.S., Construction of Orthonormal Wavelets, in "Wavelets: Mathematics 
and Applications", J.J. Benedetto and M.W. Frazier Eds., CRC Press, Boca Raton 

[22] Yoshioka D., Ground state of the two-dimensional charged particles in a strong mag- 
netic field and the fractional quantum Hall effect, Phys. Rev. B 29 (1984) 6833 

[23] Yoshioka D., Lee P. A., Ground state energy of a two-dimensional cherge- density- wave 
state in a strong magnetic field, Phys. Rev. B 27 (1983) 4986 

[24] Zak J., Dynamics of Electrons in Solids in External Fields, Phys. Rev., 168 (1968) 
686 

[25] Zak J., The kq-representation in the dynamics of electrons in solids, Solid State 
Physics, H. Ehrenreich, F. Seitz, D. Turnbull Eds., Academic, New York (1972), Vol. 
27 

[26] Zhu X., Louie S.G., Wigner crystallization in the fractional quantum Hall regime: a 
variational quantum Monte Carlo study, Phys. Rev. Lett., 70, (1993), 335 



22 



